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ABSTRACT 

We construct a nonlinear kinetic equation and prove that it is well-adapted to describ 
general multidimensional scalar conservation laws. In particular we prove that it is well-poset 
uniformly in e - the microscopic scale. We also show that the proposed kinetic equation it 
equipped with a family of kinetic entropy functions - analogous to Boltzmann’s microscopic 
77-function, such that they recover Krushkov-type entropy inequality on the macroscopic 
scale. Finally, we prove by both - BV compactness arguments in the multidimensional case 
and by compensated compactness arguments in the one-dimensional case, that the local 
density of kinetic particles admits a ‘continuum’ limit, as it converges strongly with e j 0 to 
the unique entropy solution of the corresponding conservation law. 
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1. INTRODUCTION 


Consider the scalar multi-dimensional conservation law 

(1.1) + 0)] = o, o MWSx«+, 

with given initial conditions l i(x,f = 0) = u 0 (x). We are concerned here with a Boltzmann- 
like kinetic equation which describes (1.1), as its microscopic scale, e > 0, tends to zero. 
To this end we introduce a scalar function, / e (x, v, i), which can be viewed as a microscopic 
description for the density of particles located at (x,i)ei2f X Rf with speed veR . Starting 
with given initial distribution, / e (x,u, 0), our kinetic model evolves according to 

(1.2a) [d t + a(v) • d s ]f e (x,v,t ) = j[x«.(*,t)(t>) “ /e( x > u » 01- 

Equation (1.2a) tells us that the particles are transported along 

a(u) • d x = 53 a i{ v )jT ~ i a *( ) = 4()> 

•=i ax » 

and that their collisions are governed by the nonlinear kernel on the right. Here, 

(1.2i>) u e (x ,t)= f f c (x,v,t)dv, 

J V 

denotes the local density of particles at a given (x, t) location, and the < equilibrium function , 
X«.(x.t)(v), is the signature of u e (x,t), i.e., 

, v / \ f sgnu, if (u — v)v > 0, 

( L2c ) “ { 0, if (u - v)v < 0 - 

The classical example of a kinetic model is of course the Boltzmann equation [1]. Equation 

(1.2) is closely related to the B.G.K. model of Boltzmann equation. Existence theory for 
Boltzmann equation and its simplified B.G.K. model can be found in [6], [11], respectively. 
In both cases, however, the question of convergence of the macroscopic moments to weak 
solutions of compressible Euler equations is still an open problem. (Consult [3] regarding 
an affirmative answer to this convergence question in the case of strong solutions.) In this 
paper we restrict our attention to the simpler scalar case, and we show that the proposed 
kinetic equation (1.2) is well adapted to describe strong as well as weak solutions of (1.1) as 
e 1 0. 

The paper is organized as follows. In Section 2 we show that the kinetic equation (1.2) is 
well-posed in L°°(Rt ; L r (R d x x R v )). Next, we borrow our terminology from the framework 
of Boltzmann’s kinetic equation. The microscopic scale, e, in (1.2) can be viewed as the 
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mean free path. In Section 3 we prove that the continuum or ‘fluid’ limit of the local density 
of particles, lim^o w*(z, f) is the unique entropy solution of (1.1). A kinetic construction of 
conservative solutions was carried out by Giga and Miyakawa [7]. In fact their construction 
is nothing but a fractional splitting solution of our kinetic equation (1.2), namely, a kinetic 
approximation is constructed by a succession of small time steps, in which we first transport 
and then project the particles distribution according to (1.2). Here we improve on [7] by 
identifying the underlying kinetic equation which corresponds to (1.1). It is also shown here 
that this kinetic equation is equipped with (a family of) kinetic entropy functions which play 
an analogous role to Boltzmann’s H -function. In particular, Krushkov entropy inequality 
[8], [9] is recovered in the ‘fluid’ limit £ J, 0. 

In Section 4 we revisit the question of the ‘fluid’ limit in the case of one-dimensional 
kinetic model. Here we show that the compensated compactness theory of Tartar [13] can 
be adapted as an alternative approach for providing an affirmative answer to the question 
of macroscopic convergence. The compensated compactness arguments allow us to pass 
to the continuum limit with minimal L l f] L°° information about the distribution function, 
/*, which may still oscillate around the ‘equilibrium function’ Xu- Finally, in Section 5, we 
indicate the extension of our results to the inhomogeneous case, in the presence of a (possibly 
stiff) source term. 

2 . THE KINETIC EQUATION IS WELL-POSED 

Let us rewrite (1.2) in the form 

C 2 - 1 ) g- t fe + a(v) • —f. + -f e = -Xu. , 

separating between its linear part on the left and its nonlinear kernel on the right. By 
Duhammel’s principle, (2.1) admits the following equivalent integral representation 

(2.2) f e (x, V , t) = e"« / e ( X - ta(v), V, 0) + 1 J e (T_£)/e Xu.(*-(t-r)a(v),r)(u)^T. 

The question of existence of a kinetic solution of (1.2) is now transformed into that of a fixed 
point solution for the RHS of (2.2). Fixing T,T > 0, we seek such a fixed point solution 
in 2/ oo ([0,7’]; f/ 1 (R* x Rv)). To this end, we let f e and g e be two different solutions of 
(1.2a) with corresponding densities u e (x,t) = J v f e (x,v,t)dv and w e (x,t) = f v g e (x,v,t)dv. 
By (2.2), their difference does not exceed 

\\fe(x,v,t) - g e (x,v,t)\\ LHliixRv) < e~ t/e \\f e (x,v,0) - g e (x,v } 0)\\ L i {R i xRv) + 

H / ^ ||Xu,(x,r)(' u ) — Xui,(x 1 t)(' u )||.£, 1 (.R2x.R.) < ^ 7 " - 

S Jr - 0 
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Using the properties of the signature function, %, we therefore conclude 
\\fe( x , v, t ) — g e (x,v, ~ e ^ e \\ /«(*> v >0) — g e (x, 

(2.3) 

+(! - e t/e ) $™\\M x ,v,t) - g e (x,v,T)\\ LliRixRv) . 

The inequahty (2.3) shows that the fixed point iterations 


(2.4) / e m+1 (s, V, t) = e-‘/ e / e (® - te(v), V, 0) + - e^- t )/‘ Xtir(x _ (t _ r)a(t))>r) (t,)dr > 

£ *'T=0 


are contracted (with a contraction factor of 1 — e */ e ) to a fixed-solution solution of (2.2). 
Moreover, by (2.3) this kinetic solution is unique and continuously dependent on the initial- 
data, for 


(2.5) ^max^ ||/ e (x, u, t) 0lln 1 (J*2x.R») — H/efaj^jO) 9c ( x > u > 0)||L l (flgxiiv). 

We summarize this by stating 

THEOREM 2.1. The kinetic model ( 1.2) is well-posed in L°°(Rf ; L X (R^. x /?„)). Moreover, 
its solution operator is nonexpansive in this topology, i.e., (2.5) holds. 

We close this section with several remarks. 


1. L°° - bound . To see that the solution operator associated with the kinetic model (1.2) is 
uniformly bounded, we use (2.2), obtaining 


( 2 . 6 ) 


\\fc(‘> V >t)\\L°°(Ri) < e ^ e ||/e(‘, V, 0)||£oo( R d) + 

+ (1 - e~ t/e ) max IIXu^.-t)^)!^**) < 


< e ^ e ||/ e (-, v ,0)||i / oo(flj) + 1 — e e . 


2. Finite speed of propagation . We assume that initially, f e (x, •, 0) has a compact support in 
Rv- Let us first show that / e (x, •, t) remains compactly supported. Indeed, by (2.6), /*(•, v, t) 
and hence u e (-,t) are uniformly bounded, and therefore the contributions of Xu(-,-)( v ) on the 
RHS of (2.2) are supported by ue[— «<»], where = ||u e (x, i)llx,“(j^xR^)- Consequently, 
f e (x,-,t) given in (2.2) remains compactly supported for all t > 0, with support contained 
in supp„/ e (x, •, 0) U[— Uoo, U qo]. (Note that after an initial kinetic layer of order 0(e ) , the 
contribution of the initial data in (2.2) decays exponentially fast. Thereafter, / e (x, -,t) is in 
fact ‘essentially’ supported in [— Uoo, Uqo]). 
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With this in mind we now turn to prove the finite spatial speed of propagation. We shall 
need a refined version of estimate (2.6). To this end we first observe that according to (2.2) 
- which we rewrite as 

f e (x,v,t ) = e~ t/e f(x - ia(u), v, 0) + (1 - e~ t/e ) 

f e (x,v,t) is given by a convex combination of f e (x — ta(v),v, 0) and x« r)o(v),r)(v). By 
Jensen’s inequality, therefore, we have for any convex function, U(f), 

(2.7) U(f e (x,v,t )) < e~ t/e U(f e (x - ta(v),v,0)) + \ f e (T_t)/e C/'(xu.(x-(t-r)a(«),r)(u))^r. 

c Jt — 0 

In particular, consider the case U(f ) = |/| p . If we let a „ denote the maximal speed of 
propagation, 

(2.8) aoo = {max|oi(u)| , v e supp „/ e (x, •, t)}, 
and B[r] = [— r, r] d C R d , is the ball of radius r, then (2.7) implies 

l|/e(‘> U »*)llz,»(B[r]) ^ e-t/e ||/=(') V » 0 )llL<»(B[r+ta oo ])+ 

+ “ / e(r_<)/e |IXu.(-.T)(w)||L~(B[r+(t-T)a 0O ])rf'r • 

£ Jt = 0 

Integrating the last inequality w.r.t. v we find 

ll/e(') V > i )llL«’(B[r])^ U — e J v ll/ e (’> V >0)|li=>(B[r+ta oo ])^ U + 

(2.9) 

+ (1 - e~ t/c ) max jf ||/«(-, u, r)\\ L ^ [B[rHt _ r)aoa]) dv. 

If, on the one hand, we take the p-root of both sides and let p f oo, we obtain 

(2.10) ||/ e (x,'M)||L~(B[r]xJM < maxdl/^I.U.OjIlLoofBfr+tacoJxB.),!}, 

in agreement with what we had before, consult (2.6). If on the other hand, we set p = 1 in 
(2.9), we find that the function F(t), 

F(t ) = f ||/ s (-,u,r)||Lo 0(B[r+(t _ T ) 0oo] )du, 0 < t < t, 

J v 

satisfies 

m < e- ,/ ‘F(0) + (1 F( t), 

and hence F(t ) < F( 0). Thus, we have a finite speed of propagation (< a^) of the uniform 
bound on the moments 

(2.11) [ ||/ e (x,u,t)|| L oo( B [ r ])du < / ||/ e (x,u,0)||i,oo( fl[r+tOoo]) dt;. 

Jv J v 


St= O e ^ ^ Xut(x-(t-T)a(v},T){ v )dT 

/t=o el T_t ) /e dr 
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In particular, the local density is uniformly bounded by the initial data, 

||u e (x, Olli,<»(ilgxfl+) — / ll/eC^j^jOJIlx, oo (R d)du. 

J v 

In summary, we conclude that for initial data f e (x , v, 0 )eL 1 (R v ] L°°(R f )) which are compactly 
supported in R v , the corresponding kinetic solution f e (x,v,t) remains compactly supported 
in R„ and is uniformly bounded in L 1 (i2 v ; due to a finite speed of propagation 

< Ooo, given by 

a TO = {rnax|a t (u)| I M < ||/«(x > V > 0)|| £ l( J | w;J ;»( flS )) , V e supp„/ e (s,u,0)J . 

3. Monotonicity. The signature function % u (u) is an increasing function of u. Consequently, 
the fixed-point iterations (2.4) show that 

(2.13) f e (x,v,0) > g e (x, u,0) =>■ f e (x,v,t) > g e (x,v,t), for all (a:, v), 

namely, the solution operator associated with the kinetic equation (1.2) is monotone . In 
particular, if we compare a given kinetic solution (compactly supported in R * x Ry) with 
the steady state solutions Xcon»t.( u )> i.e., if initially we have 

( 2 - 14a ) Xk{v) < f e (x,v, 0) or f e (x,v,0) <xk(v), 

we obtain 

( 214fc ) Xk(v) < f e (x,v,t) or f e (x,v,t) < Xk(v), 

in agreement with (2,6). And, since the kinetic solution operator is also conservative, the 
Crandall-Tartar lemma [5] implies the j^-contraction stated in (2.5). In fact, at this point 
we can state a little more, namely, 

4- Z/ 1 - contraction revisited . Taking into account the finite-speed of propagation, we can 
repeat - along the lines of Remark 2, a localized version of estimate (2.3) which sharpens 
the L 1 -contraction estimate (2.5) into 

(2.15) II f e (x,V,t) - 0c(z,V,t)||£*(B[r]xA.) < ||/ e (®, V, 0) - g e (x, V, 0)||Li(B[r+ta oo ]x«.)- 

5. The various estimates quoted above indicate that after an initial layer of order 0(e), the 
kinetic solution asymptotes to the ‘equilibrium function’, Xu (*,t)( w )> w h ere — as will be shown 
in the next section, u(x, t ) is the unique entropy solution of (1.1). 
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3. KINETIC ENTROPY FUNCTIONS 

Our analysis of the kinetic model (1.2) hinges on the construction of certain kinetic 
entropy functions. A kinetic entropy function in this context is a function, H(f), such that 
as in Boltzmann’s -Theorem, any solution of (1.2) obeys the additional entropy inequality 

(3.1) j{dt + a(v) ■ d x ]H(f c )dv < 0. 

We shall construct a family of such kinetic entropy functions depending on extra fixed 
parameter k,k real. To this end, we integrate (1.2a) against sgn(/ e ’X.k) over the phase 
space. Invoking a standard regularization argument of the signum function we obtain 

(3.2) j[d t + a(v) ■ d x ]\f e - Xk\dv = - i f sgn .(/, - »)(/« “ Xu.)dv. 

Noting that the expression on the right is upper-bounded by 

j sgn(/ e - X k)(fc - Xn.)dv = - Xfc| + sgn(/ e - Xk){Xk - Xu.)]dv < 

< --[/ |/e - Xk\dv ~ |u, - fc|], 

£ Jv 

we arrive at 

THEOREM 3.1. For any solution f t eL°°{Rt\L 1 {R d x x R v )) of the kinetic model (1.2), the 
following inequality holds 

(3.3) J[d t + a(v) ■ d x ]\ f e - Xk\dv < -~[J^ \ f, - Xk\dv - \u e - *|]. 

Now, the RHS of (3.3) is clearly nonpositive for 

(3.4) Xk ^ dv - L^ e ~ = _ e* Ue ” ^ ' 

Consequently, Theorem 3.1 yields 

COROLLARY 3.2. For any k,k real, the following functions 

(3.5) H(f e ) = H k (f e ) = |/ e - Xk\ 
are kinetic entropy functions, i.e., we have 

(3.6) J[d t + a(v) ■ d x ]\f e - Xk\dv <0. 
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Let us point out that our kinetic entropy functions, #*(/*), are intimately related to the 
entropy functions used by Krushkov in [8]. Indeed, as e j 0 we expect (and later on prove) 
that / e approaches With this in mind, the inequality (3.6) turns into Krushkov’s entropy 
inequality [8] 

(3.7) ii\ u - +]£ [sgn(u - k)(Ai(u) - Aj(fc))] < 0, for any real A:. 

To make this last point more precise, we shall need several lemmata. We start with 

LEMMA 3.3. Let f e eL°°(Rt ] L 1 (R V ] L°°(R d ))) be the solution of the kinetic equation (1.2), 
subject to given initial data / e (x,u, 0) which are compactly supported in R v . Assume that 
u e (x,t ) = /„ fe(x,v,t)dv satisfies 

(3.8) u e (s,0) — ► u 0 (x) in L 1 (R £), 
and 

(3.9) a subsequence of u e (x,t ) — » -u(x, t) in ££W;£‘«))- 

Then the sequence u e (x,t) converges strongly in L^ c (Rt\ L 1 (R^ i )) to u(x,t), which is the 
unique entropy solution of the conservation law (1.1), i.e., (3.7) holds. 

Note : If we take k > ||ti||jr i oo( fl d x [ 0 ,T]) ) then the entropy inequality (3.7) yields 

Wt u + 1 ~ 

i.e., u(x,t ) is a supersolution of (1.1); similarly, taking k < — ||u||£,oo(^ x [ 0i r]) shows that 
u (x, i) is a subsolution of (1.1). Hence, (3.7) implies that it(x,i) solves the conservation law 
(1.1). 

To prove Lemma 3.3 we first prepare 

LEMMA 3.4. Let f e eL°°(R(] L 1 (R V ; L°°(R d ))) be the solution of the kinetic equation (1.2), 
subject to given initial data / e (x,u, 0) which are compactly supported in R v . Then for any 
k,k real, we have 

(3.10) f | f e - Xk\dv - \u e - k\ — »0 in L^Rf x R d x ), 

Jv ej.0 

(3.11) f a(v)\f e - Xk\dv - sgn(u e - k) [ a(v)(f e - Xk)dv -» 0 in Lj^fRf x R d ). 

Jv Jv ej.0 


7 


PROOF. The vanishing limit in (3.10) follows from the inequality (3.3), for 

( 312 > T - 

0< [ f [f \fc- Xk\dv ~ \ u e - k\]dxdt < -e [ f f[dt + a(v)d T ]\f t - Xk\dvdxdt — ' 0. 
JO Jx J v Jo JxJv c| 0 

To prove (3.11) we write 

/ |/ e - Xk\dv - |u e -k\ = 

Jv 

(3.13) - f sgn (/* - Xk)(fc ~ Xk)dv - sgn(u e - k) f (/ e - Xk)dv = f (/ e - Xk)s(v)dv. 

J v J v J V 

Here, s(i>) = s(u; x, t ) is the characteristic function given by 

s(v) = sgn (f e (x, v, t ) - Xk(v)) - sgn (u e (x, t ) - k). 

Now, since s(u) is supported on the set 

V = {u|sgn(/ e - Xfc) ± Sgn(u e - A:)}, 

and since 

sgn (f e - Xk) ‘ s(u) = 2, for veV, 
we can rewrite (3.13) in the following form, 

(3.14) [ |/ e - Xk\dv - \u e - k\ = [ I fe - X*|sgn(/ e - Xk)s(v)dv = 2 t |/e - Xk\dv. 

In view of (3.10), the identity (3.14) implies 

(3.15) [ \f e -Xk\dv—>0 in L} x (Ri x R d x ). 

JvcV ffj.0 

We conclude by noting that 

/ a ( v )\fc - Xk\dv - sgn(u e - k) / a(v)(f e - Xk)dv = 

J v Jv 

(3.16) = / a(u)(/ e - Xk)s{ v ) dv — 2 / a(v)\f e - Xk\dv < 2a oa [ \f e - Xk\dv, 

Jv JvcV JvcV 

and (3.11) follows from (3.16) together with (3.15). . _ D 

Equipped with Lemma 3.4 we turn to the 

PROOF (of Lemma 3.3). By our assumption (3.8), there is a strongly convergent subse- 
quence (still denoted by) u e (x,t ) — » u(x,t). Utilizing (3.10) we obtain 

I/. - = |«» - *1 = l u - *1- 


( 3 . 17 ) 



Here the overbar denotes the weak *L°° -limit of the indicated quantities after extraction of 
appropriate subsequences, if necessary. (We note that the existence of the weak *L°° limits 
here and below are justified, since in view of (2.11), f e (x,v,t) remains compactly supported 
in R v and uniformly bounded w.r.t. e in L}(R V \ L°°(iZ^))). 

By (1.2) we have 

fe ~ Xu. = ~e[d t -I- a(v) ■ d x ]f e — * 0 in V, 

e|0 

and hence by (3.9) 

(3.18) 7 ~a~(y)f^dv = [ a,i(v)xu.dv = Ai(u e ) = A^u). 

J v J V 

This together with (3.11) gives 

/ ai{v)\fe - Xk\dv = sgn (u e - k) j a;(t/)(/ e - Xk)dv = 

J X) J V 

( 3 ' 19 ) = sgn(u - k) • ( / ai(v)f c dv - t a,(v)xfcdu) 

J V J V 

= sgn(u - fc)(A i (u) - Ai(A:)). 

Hence, in view of (3.17) and (3.19), the weak limit of (3.6) recovers the entropy inequality 

(3.20) *|u - *| + £ ^-[sgn(u - fc)(Ai(u) - Ai(*))] < 0. 

The above argument shows that the strong limit of any subsequence of u e satisfies the 
entropy inequality (3.20). Since the entropy solution of (1.1) subject to initial conditions 
(3.8) is unique, we conclude that lim c 4 ou e (x,f) = u(x, £) as asserted. □ 

We now turn to show that the continuum ‘fluid’ limit of the kinetic equation (1.2) exists 
and is governed by the conservation law (1.1). By Lemma 3.3 it remains to show that 
u e (xyt) is precompact in L 1 (i2f)). In this context there is (by now) a standard 

procedure, e.g., [4], which is based on uniform Bounded Variation (BY) estimate for each 
fixed t y coupled with equicontinuity (typically, Lipschitz continuity), in time. This brings us 
to our next lemma which states 

LEMMA 3.5. Assume that 

\\f'{x,v,0)\\ mRixLl{Rv)) = sup [- d-r / [ \f e (x + Ax, v,0) - f e (x,v,0)\dvdx] 

|A *|*0 \Ax\ JxJv 

is bounded uniformly in £. Then the corresponding kinetic solution, f e (x,v,t), satisfies 

(3.21) \\fe{x,V,t)\\ BV (BixlA(Ry)) ^ ll/e^i^jOjIlBv^xL 1 ^.))- 
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Moreover, if f e (x,v,0)eL 1 (R v \ BV(R^)) are compactly supported in R^, then we also have 
for fi,< 2 > 0, 

(3.22) ||u e (l,t 2 ) — «e(*>il)||tl(R<) < \h — ill • Goo • |l/e( x > v > 0) ||fi V(il 2 xZ, l (fl«)) • 


PROOF. Since the kinetic model (1.2) is translation invariant in spatial variables, we can 
apply the L 1 -contraction (2.5) to f e (x,v,t) with g e = f e (x + Ax, v, t) and obtain (3.21). 

Integration of the kinetic equation (1.2) over the phase space yields 

d , .v . A d 


— u e (x, 0 + £ ^7 / ai(v)f e (x, v, t)dv = 0, 


and since f e (x,v y 0)eL 1 (R v \ BV{R*)) C L 1 (i2 v ; L°°(i2f)) is compactly supported in R vy we 
may use the finite speed of propagation bound in Section 2 to conclude 

rt 2 ^ Q p 

|K(®, h) -“«(*) il)llil(i^) = / || XI -R- / a i( V )fe( x > v > T )\\L'(Ri)dT < 

J T—t\ ^ ^ vICj J V 

r f 2 

— I a oo\\fc{ x > v i' J ~^\\BV(R*xL 1 (R v ))dT' 

J r=ti 

Also, since f e (x,v,0)eL l (R v ; BV(R £)) C BV(R f x L 1 (i2 v )), the last inequality together with 
(3.21) imply the Lipschitz continuity in time, (3.22), which completes the proof. □ 

Remark . In the course of proving Lemma 3.3, consult (3.18), we established only the weak 
*L°° convergence of the spatial fluxes. However, equipped with the BV setup of Lemma 3.5 
we are able to derive strong convergence. Indeed, one may utilize the integral representation 
(2.2) to conclude that in this case we have 

f e (x,v,t) -► Xu (x,t)(v) strongly in ^([O.T] x R d x x Rf). 

This together with the finite speed of propagation imply 

II / o.i{v)f e {x,v,t)~ A(u(M))lU l ([o.rixR 2 ) < 

Jv 

^ a oo * ||/e(®jV>0 ~ Xu(x,t)('^)|| L 1 ([O l T]x ^ 1 

in contrast to the weak convergence stated in (3.18). We shall omit the details (consult 
Theorem 3.7 below), and we turn now to summarize our results by stating the following. 

THEOREM 3.6. Suppose f e (x y v J 0)eL 1 (R v ]L 1 f)L°°(R*)) such that 

(3.23) u c (x,0) = f f e (x y v y 0)dv -* u 0 (x) in L l {R*). 

Jv 
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Then the local density of the corresponding kinetic solution, u e = J v f e (x } v, t)dv, converges 
to the unique entropy solution of (LI), i.e., we have 

(3.24) f f e (x,v,t)dv u(x,t) in L°°([0, T]\ 

Jv 

and the entropy inequality (3.7) holds. 

PROOF. We begin by first assuming that f e (x,v, 0) is compactly supported in 
BV(R*)), uniformly w.r.t. e. By Theorem 2.1 (consult (2.12)), u e (x,t ) are uniformly 
bounded, and by (3.21) they have uniformly bounded spatial variation, i.e., 

||lt e (l, t)H^y(^jd) < \\fe(3'>‘V)ty\\BV(R*xL 1 (R,)) — Const. 

Hence {u e (x,t), 0 < t < T} is a bounded set in L 1 f\BV(R^) and by Helly’s theorem it 
is therefore precompact in Lj oc (R £). By (3.22), ||u e (x,t)|| i i( R d) is Lipschitz continuous in 
time, and by Cantor diagonalization process of passing to further subsequence if necessary, 
(3.24) follows. By Lemma 3.3 this completes the convergence proof for compactly supported 
BV initial data. The general case is justified by standard cutoff and BV- regularization of 
arbitrary L 1 f] L°°(R^.) initial data, consult [4]. □ 

We continue with a couple of remarks. 

1. The kinetic initial layer. We observe that Lemma 3.5 supplies us with an e-uniform bound 
on the spatial variation on the microscopic scale, (3.21). The temporal variation (Lipschitz 
continuity), however, is uniformly bounded only on the macroscopic scale, (3.22). In general, 
one cannot control the temporal variation in the microscopic scale (uniformly in e), unless 
we can prevent the possibility of a kinetic initial layer in (1.2). To this end we proceed by 

2. Preparing the kinetic initial data. In order to avoid a kinetic initial layer, we have to 
bound | jf e uniformly in e and time, in particular at t = 0. Taking into account the uni- 
form bound (in e and t ) of the spatial variation, (3.21), it remains to bound the nonlinear 
‘interaction’ kernel on the right of (1.2), \{xu. - /«)• In particular, we therefore need 

(3-25) ||/«(z,w,0) - X u .(x,o)(u)IUi(j^xii.)-^ o 0 - 

Since by our assumption (3.23) we already have that 

||Xu.(x,0)(v) - X«o(x)(w)||l»(JISxJI,) = IK(z) - wo(®)IUi(/^)-^ > o °, 

the requirement (3.25) boils down to 

(3.26) ||/«(*,«, 0) - XuoIxjHIIli^xHo) - ^ 0 - 


11 



Thus, given the initial conditions u(x,t — 0) = Uq(x), we have to prepare the kinetic initial 
data, f e {x,v, 0), such that (3.26) holds. If we prepare the kinetic initial data in such a 
manner, then we can derive explicit bounds ( uniform in time) on the error ^between the 
kinetic solution and the exact entropy solution, as told by 

THEOREM 3.7. (Error bound). We consider kinetic inital data, f e (x , u, 0)ei/ 1 (i2 t; ; BV(R*)) 
which are compactly supported in Ru. Suppose we prepare the kinetic initial data so that 

(3-27) ||/ e (x,v,0) - X«o(*)( v )IU‘(fl2xH.)-^ o °- 

Then the following error bound holds 

II fc(x,v,t) - Xu.(x,e)(w)IU‘(K-xfl.) < 

(3.28) < 2ea oo ||/ e (x,T7,0)|| il(Hv . B v(i^)) + 2||/ e (x,v J 0)-Xuo(x)( v )llLi(RixR.)-^ o 0 - 

Consequently, we have 

(3.29) / e (x, v, t ) -► Xu(x,t)( v ) strongly in L°°(Rt ; L^R? x R „)). 


Note . Preparing the kinetic initial data according to (3.27) is a strengthened version of our 
assumption (3.23). In this case, the kinetic distribution converges strongly and uniformly 
in time, to the equilibrium state Xu> as expected. Also, all the weak limits indicated in 
the proofs of Lemma 3.3 and 3.4 are in fact strong ones; in particular we now have strong 
convergence of the corresponding fluxes 


[ a t (v)f e dv ^(u) in I°°(fl+; L\R* x R v )), 

Jv 

compared with (3.18). 

PROOF. Since the kinetic model (1.2) is translation invariant in time, we can apply the 
L 1 -contraction (2.5) to f e (x,v,t) with g e = f e (x,v,t + At) and obtain 

d d 

(3-30) II xH v ) ^ || — 0)||l 1 (^xH,)* 

The kinetic equation (1.2a) enables us to upper bound the RHS of (3.30), namely, 


(3.31) 


'dt 


f e (x,v,t = 0)|| L i (J j 2xfl . ) < 


< ||[a(t/) • d x ]f e (x,v,t — 0)] L i (j^xfl.) + ~||Xu,(x,t=o) - fc(x,v,t - 0)|| i i(Hd xfl ,)- 
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The first term on the right of (3.31) does not exceed 

||[a(u) • d x ]f e (x,v,t = 0)|| L i (J ^ xJi . ) < aoo||/e(«,v,0)|| i i ( ^. BV( ^ )) ; 

the second term is less than 

1 2 

-||Xu.(x,t=o)(v) - f e (x,v,t = 0)|| L i ( ^ xiM < -||/ e (x,v,0) - XuoWIUm^xH.)- 

t c 

Substituting the last two estimates into (3.31) we end up with 

Q 

(3.32) < eo oo ||/ e (x ) u,0)|| L i (Hv . fly(B d )) + 2||/ e (x,u,0) - X«o(*)( v )l \v(RixR v )- 

Finally, we use the kinetic equation (1.2a) once more, obtaining 

II Mx,v,t) - Xu.(*,t)(v)IU l (*2x*.) < 

( 3 - 33 ) < e||^/ e (a:,u,0|| B i(ii-xR.) + e||[ a ( u )-^x]/ e (a;,v,i)|| L i (BixRi() < 

< 2EOOO ||/ e (x, V, + 2||/ e (x, V, 0) - X« 0 (*)( v ) II XA(uj *R*)> 

and (3.28) follows. 

By Theorem 3.6 we also have that u e — u and consequently that Xu 9 ~ Xu converges 
strongly and uniformly in time to zero, and by adding this to (3.28) we obtain (3.29) as 
asserted. □ 

We note in passing that the L 1 -contraction and the related BV estimates stated in Section 
2 and Lemma 3.5 are not identical with the usual L 1 -contraction statements concerning 
viscosity regularizations of entropy solutions of (1.1). In fact, at any fixed time level, we 
have 

Il/ e - ^IIlM^xR.) > / I / (/« - 9e)dv\ = ||u e - W e \\ L i( Ri) . 

By (3.29), however, the two statements coincide in the limit as e J. 0, 

II/- - 9c\\l\ JligxR,) I IXu- Xw\ dvdx = ||u - H| L i(i£)> 

and we recover the Z 1 -contraction (and the corresponding BV estimates) for entropy solutions 
of the conservation law (1.1). 

We close this section by calling attention to a rather unusual result in the theory kinetic 
equations. Namely, if u(x,t ) is a smooth solution of the conservation law (1.1), then the 
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equilibrium function Xu (*,t)( u ) an -^-solution of the corresponding kinetic equation (1.2). 

That is 


THEOREM 3.8 (exact solutions). If u(x,t)eC ^([OyT] x Rf) satisfies the conservation 
law (1-1), then Xu(x,t)(' u ) is a kinetic solution (1.2) on R* x [o,n 

Note . Theorem 3.8 is no longer valid when u(x,i) contains shock discontinuities. After the 
formation of shock waves, the corresponding kinetic solution has a “multivalued” form, e.g., 
Xu l (x.t)(w) + X[ua(x,t),u s (x,0]( v )> as m the transport collapse method of Brenier [2]. 


PROOF. We have to show that f e (x,v,t) = Xu(x,t)( u ) satisfies the kinetic equation (1.2a), 
i.e., that for any C” test function <f>(x,v, t ) 

(3.34) [ f Xu (x,t)H[dt + a ( v ) ■ d x )(f>{x, v, t)dxdvdt - 0. 

JO JRixRv 

Since the integration in R v is compactly supported (on [— Uoo> u oo])> it is enough to consider 
successively f>(x, v , t ) = <f>(x, t ) • {1, v, u 2 , . . .}, in which case (3.34) amounts to the equivalent 
conservation laws 


(3.35) 


d_ 

dt 


w 


\x,t) 


d 3 

+ s£i. 


vi p ~ 1 ai{y)dv 


= 0, p=l,2,... in T>' . 


Indeed, (3.35) are the usual entropy equalities satisfied by continuous solutions of (1.1), 
but violating (for p > 1) the Rankine-Hugoniot conditions after the formation of shock 
discontinuities. 


4. MICROSCOPIC OSCILLATIONS AND COMPENSATED COMPACTNESS 


(4.1) 


In this section we deal with the one-dimensional scalar conservation law 

du dA(u) 


dt + 


dx 


= 0 . 


The corresponding underlying kinetic equation reads 

(4.2) + a(v)£]/ e (x, v, t ) = j[xu.(x,t)(v) - f*(x, v, <)], a(-) = A'( ), 


and we raise the question of convergence of the local ‘particles density’, u e (x, t ) = J v f e (x, v, t) 
dvy towards the entropy solution, u(x , <), of (4.1). In this section we give an affirmative answer 
to this question, which is independent of compactness arguments, i.e., the BV estimates used 
in Lemma 3.4. Instead, we appeal to compensated compactness arguments, specifically, we 
employ Tartar’s div-curl lemma [13]. In this context, it is instructive to see how oscillations 
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which persist on the microscopic scale are ‘compensated’ in a manner which enables us to 
pass to the limit on the macroscopic scale. We have 

THEOREM 4.1. Let f e eL°°(Rf ; L'(R V \ L 1 f| L°°(R X ))) be the solution of the kinetic equa- 
tion (4-2). Then u e (x, t ) = f v f e (x , v, t)dv converges strongly in L p loc (R x x Rf),p < oo, to the 
unique entropy solution of the nonlinear conservation law (4- 1)- 

Remark. The conservation law (4.1) is nonlinear in the sense that there exists no interval 
on which the flux A(u) is linear, i.e., A'{u ) ^ Const. 


PROOF. Integration of (4.2) over the phase space yields 
(4.3a) d t u e + dxj v a(v)f e dv = 0. 


The corresponding entropy inequality reads 


(4.36) d t j |/ e - Xk\dv + d xj v a(v)\f e - Xk\dv < 0. 

Since by (2.11) the left-hand side of (4.3b) lies in Murat’s lemma [10], [13] implies 

that the negative measure on the right of (4.3b) lies in a compact set of H loc (R x x Rt )• Hence 
we can apply the div-curl lemma [13] to the left-hand sides of (4.3a) and (4.3b), which gives 

u e 1 \ J e - XJtl^ - / a(u)/ e du- / |/« - Xk\dv = 

Jv Jv 

m n 

= u e -J a(v)\f e - Xk\dv — J a(v)f e dv ■ JJ\f e - Xk\dv- 

We recall that the overbar denotes the weak *L°° -limit of the indicated quantities after 
extraction of appropriate subsequences, if necessary. Following [12], we can rewrite (4.4) in 
the equivalent form 

(4.5) ( u e -Ue)- 1 <*(«)! fe - Xk\dv = Jjfe- Xk\dv ■ a{v)f e dv - £ a(v)f e dv). 


Using (3.10) and (3.11), the last equality is further simplified into 

(4.6) ( u e — u e ) ■ sgn(u e — k) J a(v)(f e — Xk)dv = |u e — A:| - a(v)f e dv — J f a(v)f e dv). 

We now examine (4.6) at an arbitrary fixed location (x,t); with k = u e (x,t ) we find after 
little rearrangement 


(4.7) 



= 0 . 
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Of course, by (4.2) 

fe ~ Xu. = ~e[d t + a(v)d x ]f e — * 0 in V, 

* io 

hence 

( 4 - 8 ) / a(v)f e dv = I a(v)xu.dv = A(u e ); 

Jv Jv 

Also, we recall with k = u e (x, t ) we have 

( 4 - 9 ) f a(v)xkdv = A(u e ). 

J v 

Inserting (4.8) and (4.9) into (4.7) we find 

\u e - u e | • (A(u e ) - A(u e )) = 0. 

This implies that 

( 4 -10) A{u e ) = A(u e ), 

for otherwise, |u e — u e |(i,t) = 0, which in turn leads again to (4.10). Taking the weak limit 
of (4.2), we obtain with the help of (4.8) and (4.10), 

+ ^A(u.) = 0 . 

Thus, (a subsequence of) u e (x, t ) converges to a weak solution of the conservation law (4.1). 
Moreover, in view of the nonlinearity of A(u), equality (4.10) implies that u e (x,t ) converge 
strongly in L p loc {R x X R+),l < p < oo, consult Tartar [13]. Using this fact together with 
Lemma 3.3 we conclude that u e converges strongly in Lf^R,. x Rf),p < oo, to the unique 
entropy solution of (4.1), as asserted. □ 


5. CONSERVATION LAWS WITH A SOURCE TERM 


In this section we extend the above results to inhomogeneous scalar conservation laws 
(5.1) 


JjK 1 . *)! + 2 *))] = S(x, t, u), (x, t)iRi X R } , 


where S(x, t, •) is an L°°(R* x R + ; C 1 ) source term satisfying S(x, t, 0) = 0. 

The corresponding kinetic model equation reads 

(5.2) [d t + a(v ) • d x ]f e (x, v, t) = j[Xu.(*,t)( u ) ~ fe(x, v, t)] + 5'(x, t, v)f e (x, v, t), 

and is augmented with the constitutive relations (1.2b), (1.2c). 
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A uni que kinetic solution for (5.2) can be constructed, as before, by Banach fixed point 
iterations which yield 

THEOREM 5.1. The kinetic model (5.2), (1.2b-c) is well-posed in L°°(i? t + ; L\R* x R v )). 
Moreover, if f e and g e are two different inhomogeneous kinetic solutions of (5.1), and if we 
let S'^it) denote 

(5.3) *Soo(t) = | max S'(x, t,v) \ v e supp „/ e (x, v, t ) (J supp^x, v, t)| , 
then we have 

(5.4) ma < x i \\fe{x,v,t)-g e (x,v,t)\\ L i( F 4 xRv ) < e-ft=o s ^ dt -\\f e (x 1 v,0)-g e (x,v,0)\\ L ^ R d xRv y 

We shall only indicate the proof of the L°°(Rt\L l {R* x R v ) stability stated in (5.4). The 
difference between the kinetic solutions f e — g e (with corresponding local densities, u e (x,t) = 
f v fe{ x > v > i )d‘ v and w e (x, t) = /„ g e (x,v,t)dv) satisfies 

[d t + a(v) ■ d x ](f e - g.) = ^[(Xu.(,.t)(v) - *..(-.*>(")) “ (/• ~ 9c)] + S'(x,t,v)(f e - g e ). 

Multiplying this by sgn(/ e — g e ) and integrating over R v and R £ (in this order), we obtain 

^|| /.(*, t) - g e (x , u, t) \\v(RixR v ) < Slot*) * ll/.(*i 0 ~ 9*(x, ”, OIIx*(JSxr.)i 

and (5.4) follows. D 

We conclude with several remarks concerning the entropy inequality. 

The corresponding inhomogeneous kinetic entropy inequality now reads 

(5.5) J [d t + a(v) • d x ]\f e - Xk\dv < J^S'(x,t,v)\f e - Xk\dv. 

Moreover, by arguing along the lines of the stability estimate (5.4) we find that for BV(R^) 
source terms we have 

||/ e (x,V,T)|| BViRixL^Rv)) < e^=° 5oo(t) ‘ ,t ||/ e (®. u »0)l|BV(RjxLi(H.)) + 

(5.6) T 

+ / e^r=t s '^ jT • max ||5'(x, t, v)\\ B v(Ri)dt ■ ||/ e (a:,'y,0)|| L i(H|xB.)- 
Jt=o v 

This allows us to keep the convergence statement of Theorem 3.6, 

[ f e (x,v, t)dv — ^(x,*), in £°°([0, T], I 1 ^)), 

Jv ej.0 

in the inhomogeneous case (5.2). In view of (5.5), we are also able to recover the macroscopic 
‘continuum limit’ entropy inequality for the above limit u = u(x,t), which in this case 
amounts to 

l\u - k\ + X) J-[ s S n ( u ~ *0(^i( n ) “ ^(*0)1 ^ 

uZ 1=1 

< sgn(u - k)[S(x, t, u) - S(x, t , k)], for any real k. 
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